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Abstract 

We discuss the potentially important role played by preheating in certain variants of the 
curvaton mechanism in which isocurvature perturbations of a D-flat (and F-flat) direction 
become converted to curvature perturbations during reheating. We analyse the transition 
from inflation to reheating in some detail, including the dynamics of the coupled curva- 
ton and inflation fields during this transition. We discover that preheating could be an 
important source of adiabaticity where parametric resonance of the isocurvature compo- 
nents amplifies the super-horizon fluctuations by a significant amount. As an example of 
these effects we develop a particle physics motivated model which we recently introduced 
in which the D-flat direction is identified with the usual Higgs field. Our new results show 
that it is possible to achieve the correct curvature perturbations for initial values of the 
curvaton fields of order the weak scale. In this model we show that the prediction for the 
spectral index of the final curvature perturbation only depends on the mass of the curvaton 
during infiation, where consistency with current observational data requires the ratio of 
this mass to the Hubble constant to be < 0.3. 
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1 Introduction 



In a recent paper we proposed a new mechanism for generating curvature perturbations during 
reheating in hybrid inflation [1]. The idea was that, during inflation, one of the slowly rolling 
scalar fields, other than the inflaton, could develop large isocurvature perturbations, and that 
this could become converted to curvature perturbations during the process of reheating. This 
is quite different to the curvaton mechanism as originally proposed [2,3] since the scalar with 
large isocurvature perturbations is not late decaying^. Instead the mechanism in [1] relies on 
the simple observation that just after the end of inflation, at the onset of reheating, the fields of 
hybrid inflation [8,9] tend to share out the vacuum energy which dominates during inflation. This 
typically leads to hybrid flclds with energy densities of the same order of magnitude during the 
epoch of reheating, allowing efficient conversion of isocurvature to curvature perturbations during 
reheating. This mechanism should not be confused with the late decaying curvaton mechanism, 
whose physics is not related to reheating. However for ease of terminology we shall continue 
to call the hybrid field whose isocurvature perturbations during infiation become converted to 
curvature perturbations during reheating the "curvaton" , even though our mechanism is different 
from the curvaton mechanism as originally proposed. 

The role of our curvaton will be played by a D-fiat (and F-fiat) direction^ during infiation, 
coupled to the waterfall field N of hybrid infiation. Therefore it is a "coupled curvaton" , which 
again is different from the original curvaton idea. Furthermore, being a D-fiat direction, typically 
it will decay through its gauge interaction before the singlets whose decay rate is controlled by 
small Yukawa couplings. We gave an explicit example of the general mechanism in which the 
the Higgs fields of the supersymmetric standard model played the role of a coupled curvaton, i.e. 
developed isocurvature perturbations during inflation which subsequently become converted to 
curvature perturbations [1]. We discussed the nature and evolution of the Higgs perturbations 
during the epoch of inflation, and described the evolution of the perturbations during reheating 
after the Higgs decayed. After Higgs decay, the Universe consists of a mixture of matter (the 
oscillating singlet flelds) and radiation (Higgs decay products), and we were able to show that 
the isocurvature perturbations are converted into curvature perturbations in this framework. 
However in the previous paper we did not analyse the transition from inflation to reheating before 

-"^For recent discussion on several scenarios for the late decaying curvaton see Refs. [4- 7]). 

^In supersymmetric models, a flat direction has to be both D-flat and F-flat. For the inflaton and the N 
mediating field, being singlets, the D-flatness condition is trivial fulfilled. By calling the curvaton a D-flat 
direction what we want to stress is its non-singlet nature. 
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the Higgs curvaton decayed, and we also neglected the very important effects of preheating. 

The purpose of the present paper is to discuss both of the above effects in detail, with 
surprising results (at least to us). Due to the presence of the coupled curvaton, the isocurva- 
ture perturbation (entropy) between the inflaton and the N waterfall field is converted into the 
adiabatic curvature during the transition from inflation to reheating. The value of the total 
curvature perturbation at the end of the transition from inflation to reheating is then given in 
terms of the value of the coupled curvaton field rather than the inflaton. As it turns out, with 
the coupled curvaton in hybrid inflation, the initial value at horizon crossing of the isocurvature 
perturbation decreases during the transition from inflation to reheating. To be precise, we find 
that during the transition from inflation to reheating, before the Higgs curvaton decays, the 
curvature perturbation drops by several orders of magnitude from our estimate in [1]. Fortu- 
nately when preheating is taken into account the curvature perturbation may subsequently be 
increased by a similar number of orders of magnitude, resulting in an acceptable final amplitude 
of curvature perturbations, consistent with COBE and WMAP. 

The layout of the remainder of the paper is as follows. In section 2 we briefly review the 
model, and in section 3 we describe the inflationary trajectory and the transition from inflation 
to reheating. In section 4 we describe the evolution of the fluctuations during inflation, and 
we carefully analyse the effect of the transition from inflation to reheating. It is during this 
transition that the curvature perturbation flrst changes. In section 5 we discuss the effects 
of preheating on the curvature perturbation, while allowing the Higgs to decay to radiation. 
Section 6 contains a discussion of the prediction of the spectral index and its correlation with 
the parameters which are consistent with giving correct structure. Section 7 concludes the paper. 

2 The inflationary model 

The supersymmetric hybrid inflation model is based on the superpotential [10]: 

W = XNH^Ha - K(t>N^ , (1) 

where N and are singlet superflelds, H^^d are the Higgs superflelds, and A, k are dimensionless 
couplings. Other cubic terms in the superpotential are forbidding by imposing a global U{1)pq 
Peccei-Quinn symmetry. The superpotential in Eq. (1) includes a linear superpotential for 
the inflaton fleld, 0, typical of hybrid inflation, as well as the singlet N coupling to Higgs 
doublets as in the NMSSM. We notice that in standard hybrid inflation [8,9] the singlet N will 
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quickly settle to its false vacuum value at zero, leaving the inflaton flat direction to slow-roll 
until it reaches the critical value. This means that during inflation the /i term generated by 
the coupling of N to the Higgses in (1) vanishes, and the Higgses become an F-flat direction. 
Therefore, we can have an additional flat direction made by a combination of the Higgs fields 
which satisfies the D-fiatness condition [11], i.e., such that D-term contribution to the potential 
{92 +5'i)(-^u ~ -^1)^/8 vanishes; for example taking them to be equal, if„ = = h. Therefore, 
a more general inflationary trajectory is given by allowing the Higgs fields to take non- vanishing, 
although as we will see, small values. 

Inflation takes place below the SUSY breaking scale. Including the soft SUSY breaking 
masses, m^, m^r and m/j, and trilinears A,^, Ax, the potential for the real part of the fields 
reads^: 



V = V{0) + ^N' + ^h' + K\<f>-<l>t){<l>-<f>:)N' + ^N'h' + XNh'l^ 



^ ,2 i2 I ^^2i2 I ^^2 Ar2 



+ _^2^2 ^ ^ ^ ^2) 

The constant V^(0) has been added to ensure the necessary vanishing of the potential at its global 
minimum, and the critical value of 0, are given below in terms of the soft SUSY parameters: 



c 



1 ± 



2V2K \ \ Al 



4m^ 



(3) 



As we will see in the next section, during inflation (when > 0+) the inflaton fleld, 0, is slowly 
rolling along the potential and also the fields and h get small values but such that the N field 
dependent mass squared, m%, is dominated mainly by the term — 0+)(0 — 0~). Once at 
the critical value, = 0^, fh% changes sign and both singlets, and N, roll down towards the 
global minimum 0o and Nq ending inflation^. 



v/2 2k \ Al 



1-^- (5) 



^The axionic (imaginary) part of the complex fields can be set consistently to zero, and it will not affect the 
evolution of the real scalar fields either during inflation or the reheating period. 

^We shall see that h <C No,(j)o and also that rn^jirij^ ^ k^Nq. Therefore either contributions from the Higgs 
field, h, and from the bare masses, m^,7n%, can be neglected in order to calculate the global minima, (j)o and 
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The Higgs field, h, will oscillate once inflation ends around the minima which is zero as long as 
the Electroweak Symmetry is still preserving^. 

Before describing the inflationary trajectory, we briefly discuss the values of the parameters 
of the potential in Eq. (2). They are uniquely determined by embedding this model into an 
extra- dimensional framework with just one fundamental ("string") scale M* ~ lO"*^^ GeV [13]. 
Embeding all the Higgs fields and singlets and N in the bulk, while all the matter fields live 
on the brane, it is possible to show that [13,12] 

A^«^ fi^y ^10-10, (6) 
\mpj 

^«~Aa~MJ— -lO'GeV, (7) 
\mpj 

ms ~ rriN ~ ruh ~ — [ — - ) ~ 1 MeV , (8) 

An \mpj 

ViQf'^ ^mJ —\ - 10^ GeV , (9) 
\mpj 

where mp — Mp/\^ — 2.4 x 10-*^^ GeV is the reduced Planck mass. Using Eqs. (4) and (5) we 
obtain: 

00 - iVo - M, ~ lO^^GeV , (10) 
and the Hubble parameter during inflation is of the order of 

H=—^ — 10 MeV. (11) 
v3mp 

Notice that < H Ioy a — (t),N,h. 

We will use those values parameters in the following sections when presenting numerical 
results. However the mechanism is general for a potential like Eq. (2) independently of parameter 
values, and the analjd;ical estimations are presented in general in terms of the model scales^, i7, 
0c, /t0c- We will refer in the following to the "Higgs" like the coupled curvaton, D-fiat direction 
during inflation, included in Eq. (2). 

3 Evolution of the background fields. 

Including the Higgs terms in the potential Eq. (2), the term proportional to XNh"^ can induce a 
non-zero value for h and during inflation providing that the coefficient of this term is negative 

^In a future work [12] we will show how the electroweak symmetry breaking is implemented in this model 
and their cosmological consequences. However here we arc going to consider that the electroweak symmetry is 

preserved all tlic' time, since our principal c;onclusions will not be affected by this issue. 
^For order of magnitude estimations, we will take (/>c ~ (/)^ ~ (/)o ~ -/Vo, and A ~ k. 
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^ — K(f) < 0. However, in order to ensure a flat direction for the Higgs, and therefore allowing 
it to slow-roll during inflation, its effective fleld dependent mass must be smaller than H'^. The 
latter receives contributions from all the fields, 

^ X^i3e + N^) + 2XNi^-K<P)+ml<H\ (12) 

This implies that both the fields N and h must take small values during infiation, with 

N ( hV .r, h H . 

-< — ^10-^°, -< — «.10-^ 13 

(Pc \K'(PcJ (Pc >^(Pc 

SO that N <^ h <^ (f). These values are small enough not to disturb the infiaton slow-roll, 

^ K'N'+mlc^ml<H'. (14) 

The Higgs and inflaton fields evolve independently of each other and N, following the approxi- 
mate slow-roll trajectories given by: 

0(t) ~ -v^H<j){t) , h{t) ~ -VhHh{t) , (15) 

where rjh — m|/(3if^), 77^ = m^/(3if^). Inflation as such is still controlled by the inflaton 0, 
and will end when (j) reaches one of its critical values 0^. 

On the other hand, the N fleld follows the evolution equation: 

N + 3HN + Vn^O, (16) 

with 

^iv=|^ = {i^'N' + 2n'i<f>-<Pt)i<i>-<P;) + X^e)N + \h'i^-n<p) 

^ u;l(<t>)N + Xh\^~K<t>), (17) 

where in the second fine oj%{(t)) » H'^ is the fleld dependent effective squared N mass, where 
we have dropped the small term k^N"^ -\- X^h? using Eq. (13). During inflation we can take 
the slow-rolling flelds h and to be approximately constant in Eq. (17). Thus, the N fleld 
will oscillate with an amplitude damped by the expansion as oc a{t)~^/'^ and a frequency of the 
order of ojn^cI)) ~ 0{K(j)c) ^ but around a non vanishing value given by the quasi-constant 
contribution in Eq. (17), i.e. 
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Figure 1 : We show the evolution of the background fields (solid lines) and time derivaties (dashed lines) during 
inflation up to the critical value = 0^ (LHS), and from the critical value up to the global minimum (RHS). 
We have taken the values of the parameters: (/)o = 10^^ GeV, </)|j" = 2(/)o, k0o = 1 TeV, A\ — 1.5V2k(/)o, 
= rah — O-IH, and the value of the Higgs field h{0) = 247 GeV. The numerical integration starts around 60 
e-folds before the end of inflation. 

where Fosc[^n{4')] represents the oscillating part. Given that in any case the amplitude of the 
oscillations quickly decays due to the exponential expansion, after some e-folds they become 
completely negligible, and Eq. (18) gives the trajectory: 



(19) 



that is, the motion of effectively follows the valley of minima in the potential given by 
Vat = 0. This inflationary trajectory for is independent of the initial conditions, provided 
that the other fields are already slowly rolling. Although in origen A^ is a "heavy" field (mass 
larger than the Hubble rate), the dynamics is such that it behaves effectively as a "light" one. 
The quasi-exponential expansion during inflation settles it to the minimum, but the minimum 
is a time-dependent, slow- rolling one. This has implications also from the point of view of its 
quantum fluctuations: they will behave in a similar way as those of a light field. 

The evolution of the background fields and their time derivaties can be seen in Fig. (1). We 
have choosen to plot the evolution against the diference |0 — in order to detail the changes 
around the critical point, which we analyse in the remaining of this section. In the left panel, 
we have plotted the evolution during infiation up to the critical point, and in the right one that 
from the critical point to the global minimum. Time fiows in both from left to right, and we 
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have started the numerical integration around 60 e-folds before the end of inflation. Although 
the N field is also "evolving" during inflation, this is standard hybrid inflation in the sense that 
it is the N field which triggers its end as can be seen in RHS plot in Fig. (1), when (p leaves the 
slow-rolling trajectory ps; Constant. 

In the LHS plot in Fig. (1), it can be seen that close to but before reaches the critical value, 
as uj%{(f)) — > 0, the values of the derivatives of increase very rapidly and the effective "slow 
roll" of the N field given in Eq. (19) is violated. In particular, N becomes now the dominant 
term to cancel out the Higgs and quasi-constant contribution in the equation of motion, such 
that 

N + 3HN + T/jv = =^ TV ~ Xh^{K(l) - ^) . (20) 

v2 

We want to estimate the ratio (N/N) at the critical value = 0^, which will be of use later 
when computing the metric perturbations. The transition around the critical value, and from 
there to the global minimum, takes a very short interval of time with HAt -C 1, so that we 
can safely neglect the effect of the expansion in the analytical estimations. By matching the 
above trajectory Eq. (20) with the inflationary one Eq. (19), we see that the transition between 
trajectories occurs when 

and we have denoted this point by a subindex "— ". Given the short interval of time lapsed 
between this point and the critical value, practically the values of N and N at 0J are given by 
those derived from Eqs. (19) but evaluated at «;(0 — 0^)- above, which gives the ratio, 

where we have used the slow- roll approximation^ 0/0 ~ rj^H. 

''This is only consistent if by that time the value of Nc is still small enough not to disturb the others trajectories, 

i.e, 

A/i2(k0+-Aa/V2) f<j)Y^^ ( H 



" {Kc^t^ilt - 0c-))V3 ^ \n<i>t) ^^^^ 

which implies 

Xh < 0{H (^^K(t)}j )<H, (24) 

which we assume to hold hereafter unless otherwise stated. That is, end of inflation in the standard sense will 
proceed only after (f) crosses the critical value. If Eq. (24) were not fulfilled, for example when the value of the 
Higgs at horizon crossing is very close to the uppc;r limit given in Eq. (13), then inflation will terminate instead 
before (p reaches the critical value by the increasing values of N and h. 



7 



Then, as shown in the RHS plot in Fig. (1), as soon as crosses the critical value and 
the N field dependent squared mass becomes negative, the value of the field starts growing 
exponentially, 

N ~ N,exp{ f ^2/s;(0+-0)(0-0-)rft) , (25) 



until finally the large values of N modify the slow-roll of the infiaton field, and this also moves 
fast towards the minimum. Afterwards, N and will move towards a straight-line trajectory 
[14] given by 

N = v^(0+ - 0) . (26) 

Once N reaches the trajectory in Eq. (26), the slow-roll regime for the inflaton (and the Higgs) 
ends. The Higgs also gets desestabilized due to the increasing N field, and it will just become 
increasingly coupled to the others, but at a slower pace; i. e., the Higgs does not really affect 
the transition of the infiaton and the N from the critical point to the global minimum. 

4 Curvature perturbation during inflation and the sub- 
sequent phase transition. 

Our final aim is to compute the spectrum of primordial curvature perturbations, initially origi- 
nated during infiation. In the presence of several light (slow-rolling) scalar fields 0ck) the comoving 
curvature perturbation [15] can be written as [16,17] 

-^-HjziM^Q^, (27) 

in terms of the gauge-invariant scalar field amplitude perturbations, the Sasaki-Mukhanov vari- 
ables Qi [18], 

Qa = <^0a + ^V', (28) 

where 5 (pa is the quantum fiuctuaction of the background field (pa, and ip the gauge-dependent 
metric perturbation (see Appendix A). The total curvature perturbation is then given as the 
weighted sum of "individual" curvature perturbations 

TZa = hH^ , (29) 

(pa 



such that. 



'?=i:7^K„=2:^^'j., (30) 
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with = Pa + Pa for a scalar field, Pa being the energy density and the pressure. The 
RHS of Eq. (30) applies in general to a multi-component Universe, independently of the nature 
(scalar fields or other) of the components. In order to follow the evolution of the curvature 
perturbation from inflation to the reheating stage, we will follow in detail that of the individual 
components TZa and how their relative contribution to TZ changes. 

During inflation, nevertheless, the field that dominates the total curvature perturbation is the 
one with the largest velocity, this being the infiaton field (p ^ N . From Eq. (30), we simply 
have TZinf ~ TZ^. However, as we will see next, the curvature perturbations of the Higgs and N , 
TZh and TZn, are orders of magnitude larger than that of TZ^. Relative perturbations between 
different components are the origen of the so called entropy or isocurvature^ perturbations, 
TZa — T^p- Therefore in our scenario, given the strong hierarchy between the values of the fields, 
there is a large "isocurvature" perturbation between the infiaton and the Higgs, and the infiaton 
and the A^. Those can later seed the total curvature perturbation, if one of the field (kinetic) 
energy density other than that of the infiaton becomes non- negligible [19-24]. This is what 
happens at the end of inflation: as soon as the inflaton field reaches the critical value, the field 
N evolution changes and moves fast towards the global minimum, with increasing kinetic energy. 
Therefore, during the transition from the critical point 0^ towards the minimum (pQ, the held N 
with a large kinetic energy gives the dominant contribution to the total curvature perturbation 
in Eq. (30); i.e., the entropy perturbation between N and (p becomes adiabatic curvature 
perturbation. This can be seen by writting the evolution equation for the total curvature^ in 
the presence of the scalar flelds [23] 0, N and h: 

iz -3HY.{cl-4)^{na-n,) 

a,l3 "-T 



- 3Hicl-cl)^{TZ,-n,), (31) 



where as a shorthand notation we have defined: = Pa + Pa = 'Pa', and Ht is the total density 
energy, Ht = h^-\- + hh] = Pa/ Pa is the sound speed for each component, 

+ (^^' 

*So called as far as the originating field or fluid component does not dominate the total energy density, and 
therefore its fluctuations make little contribution to the curvature one. 

^As we are only interested in large scale fluctuations, we neglect the term that goes like {k/aH <^ 1) in the 
equation. 
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Clearly when N cf)^ the first term in Eq. (31) will dominate, such that TZ ~ TZ^- Notice also 



that during the transition (f) and N behave as "different fiuid" components with different sound 
speeds: the former is still slowly-rolling while N is increasing exponentially. 

For a multi-scalar system, rotating the fluctuations in field space, as done in Ref. [17], 
helps to identify all along the "adiabatic" mode responsible for the the total curvature, and 
the "entropy" modes that influence its evolution. For comparison with this approach, we have 
included the rotated equations in Appendix B. Entropy modes source the adiabatic one when 
the fields follow a curved trajectory in field space. This matches the analitical behaviour in the 
next subsection: during the transition, the trajectory is curved by the motion of N , such that 
the curvature changes. 

4.1 During inflation. 

The equation of motion for the gauge invariant quantum fluctuation Qa, with comoving wavenum- 
ber A;, is given by [25] 



inflation; i.e., their fluctuacions behave nearly as those of a massless fleld. Hence, and Qh 
will be frozen to a constant value once outside the horizon, k < aH, given approximately by the 
value at horizon crossing k — aH [26], 



with 71^ = HQ^/(f) <^ TZh = HQ^/h during inflation. 

On the other hand, neglecting for simplicity the sub-dominant^'' terms coming from metric 
contributions on the R.H.S in Eq. (33), the evolution equation for reads 

k^ 

Qn + SHQn + (^ + Vnn)Qn + Vn4>Q4, + YnhQu ^ . (35) 

Like in the case of the evolution of the background fleld N. now the large mass term Vnn ^ 
0{K'^(f)1) ^ gives rise to oscillations with an amplitude decaying as a^^, but displaced from 
zero due to the mixing with and Qh- Therefore, once the amplitude of the oscillations 
^°The metric contributions are of the order of 0{4>^/ p), and therefore neghgible during inflation. 





(34) 
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becomes negligible, the N field fluctuation will also tend to a quasi-constant value given by 

Qn^-^Q,-^Qh, (36) 

VNN VNN 

and using Eq. (19) we have for Qn and the background velocity N: 

Qn ^ 7V(^2^ + C^(0)^j c^2Ar^, (37) 

N ^ N(2^ + CMfj , (38) 

and Qh/h » Q4>/4>- Then, using Eqs. (37) and (38), TZn is given during inflation by: 

" (41) 



with 



2r//i + C^{(j))rj^ h 

which is of the same order of magnitude than TZh during inflation, while is not too close to 
the critical value. We stress that it is the Higgs fiat direction which holds the N field to a 
non- vanishing, non-decreasing value during inflation, in a kind of "slow-roll" trajectory; in a 
similar way their quantum fluctuations Qn follow Qh instead of being redshifted away. 

We do not require any special initial conditions for Eq. (40) to hold, apart from the fields 
being already in their slow-roll trajectories. Given that inflation will last more than this, the 
fields will have most probably reached that trajectory much before the largest observable scale 
today become super-Hubble, around say 60 e-folds before the end of inflation. 

Therefore, up to near the critical value, TZ^ and TZ^ are frozen to their values at horizon 
crossing, as expected for nearly massless and decoupled fields, while 7^ at in Eq. (40) decreases 
as the ratio N/N evolves in time. In order to estimate its value by the time inflation ends, we 
use the fact that the ratio Qn/N remains approximately constant, and given by Eq. (37) all 
along inflation. Therefore, the contribution TZ^ at the critical value is given by Eq. (40) but 
evaluated using {N/N)^ in Eq. (22), 

1/3 



- 2 ( ^) f ^W,. ^ 2 i-jTrJ^^-—- ) 7^,. . (42) 
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Thus, by the end of inflation the value oUZn is suppresed with respect to its value at the time 
of horizon crossing by a factor 

(^) (^) ^"^^ ' ^^^^ 

which for r]^ ~ 0.01 and — 1 TeV, with if ~ 10 MeV (0c — lO^^Gel^) would give a 
suppresion factor approximately of the order of 10~^. 

4.2 During the phase transition. 

After crossing the critical value, the ratio Qn/N still remains practically constant until N reaches 
the straight line trajectory Eq. (26). Both evolutions, that of and its fluctuation, are now 

governed by the tachyonic instability in their effective mass. This means that TZn oc Qn/N 
diminishes again by a factor N/N, before the fields reach the global minimum. Let us denote 
the point when they reach the trajectory Eq. (26) by a subindex "+". Thus, the value of 7?.jv+ 
is given by 

with Ri\fc given by Eq. (42) above, and (0 — (p^)- by Eq. (21). By using Eq. (25) and 
the condition Eq. (26), we have (0^ — 0)+ fti (0 — 4>'^)^. Therefore, during the waterfall the 
value oi TZn would only decrease at most by an order of magnitude. Numerically we found 
(0+ -0)+^ 0.2(0 -0+)_ 

In the final part of the transition from 0c to 0o, the evolution of the quantum fiuctuations 
become coupled through the interaction in the potential. Given that N — •\/2(0 — 0+), and 
N — — -\/20, for the fiuctuations we have Qn oc N and Q,f, (x (j) such that Qn — —V^Q<p [14]; 
i.e, once Eq. (26) is fulfilled, TZn and T?.^ remain constant. By the time they reach the global 
minimum we simply have: 

'7^00 = Tim = T^N+ ■ (45) 

Because the background Higgs is still moving towards the minimum but with /i <S 0, the total 
curvature perturbation is then given by the contributions of the 2 singlets, with 

7^o ^ 7^^o = Tim . (46) 
Using Eqs. (44) and (42), we then have: 
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And for the spectrum, P-ji = A;^(|7?.p)/(27r^), we obtain 

- Wi'Pt - 'Pc)'Pt ) v'J% v'JX' ^ ^ 

This can be compared to the initial value of the amplitude of the spectrum, given by that of the 
inflaton at horizon crossing, 

P^" «^ ^ , (49) 

with 

Clearly, Eq. (48) can be larger than the standard value Eq. (49), if h^jcj)^ <^ {4'cV4>I^pY^^ i 
which is required in any case to ensure slow-roll inflation. In this scenario, the total curvature 
at the end of inflation grows from its value at horizon crossing, as the flelds dynamic changes 
and the relative contribution of the N fleld to the total curvature become non-negligible during 
the phase transition. As we have already mentioned, this can be interpreted as converting the 
isocurvature perturbation between N and into the adiabatic one TZ, see Eq. (31). However, the 
presence during inflation of a large amplitude TZf^ is initially controlled by the Higgs. Because 
of that we prefer to refer to the Higgs h field as "curvaton" , instead of the mediating field N . 

The Higgs does not contribute per se to the total curvature at this point, because its velocity 
is still negligible with respect to the others. Nevertheless, the evolution of /i, h and Qh resembles 
that of the other fields. First, diiring the waterfall when N becomes larger than H'^/4>, the slow- 
roll approximation for the Higgs breaks down. The evolution of the Higgs is then dominated by 
its hnear coupling to N, such that 

h « 2XhN{K(j) - ^) , (51) 
v2 

and then h oc a/ZV. The fiuctuation Q/t follows the background field, and as a result Rh decreases. 
Until the other fields reach the the trajectory Eq. (26); then again we recover the relation oc h, 
and TZfi will tend to level with the others. 

The behavior of the spectrum of the different contributions Ra during inflation and the 
phase transition, is plotted in Fig. (2). The evolution is shown in logarithmic scale versus the 
value of — same than in Fig. (1). In the LHS plot, we follow them from horizon crossing 
k — aH, taken to be 60 e-folds before the end, through inflation, up to the critical value (strictly 
speaking, at inflnity outside the plot); and in the RHS plot we have their evolution during 
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Figure 2: Evolution of the spectrum of the different curvature perturbations for each field, versus 

for k — aH 60 e-folds before the end of inflation. On the left plot it is shown their behaviour during inflation, 

and on the right one from the critical value up to the global minimum. Values of the parameters same than in 

Fig.(l). 



the transition towards the global minimum. This plot summarises the previous findings: during 
infiation TZh and TZ^ remains constant, while the initial TZn — T^h* is suppressed as it approaches 
the critical point due to the increasing (see Fig. (1)). When infiation ends, falls fast into 
the global minimum, dragging in the way also the infiaton and to some extend the Higgs (Fig. 
(1)). As a result, their quantum fiuctuations will tend to follow each other, such that finally 
T^ho — T^m = TZ<i)0- 

Nevertheless, given the model parameter values in section 2.1, it seems that the amplitude 
of the spectrum in Eq. (48) is too small to match the COBE value, P^'^ ~ 5 x 10~^, unless 
we fine-tune the value of r/^. However, this is the value at the end of infiation. This will again 
be modified by the subsequent evolution of the fields during reheating, in particular it can be 
amplified by "preheating" effects. 



5 Curvature perturbation during reheating and preheat- 
ing. 

When the background fields oscillate around their global minimum values, all the three effective 
masses rhi become of the same order of magnitude, i.e., rhi ~ O(fi;0o) ~ 0{1 TeV), this being the 
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typical value for the frequency of the three oscillating fields. Given the mixing in the potential, 
they will end oscillating with similar amplitudes. That is, the vacuum energy term 1^(0) ends 
being equally redistributed among the three matter fields. Later on, they will decay and transfer 
that energy into radiation. This is the essence of the reheating process. In our scenario, the 
Higgs has the largest perturbative decay rate Th due to either its large top Yukawa coupling, 
or gauge couplings aw, i-e., Th — aw{i<-4>o)- On the other hand, the singlet fields are very long- 
lived, due to their tiny Yukawa couplings, k ~ A ~ 10~^°, with T^^^ ~ ^^(k^o)- Therefore, at 
least perturbatively, we expect the Higgs to be the first in being converted into radiation almost 
immediately after the oscillations start. At the same time, Higgs perturbations are converted 
into those of the radiation fluid. 

Nevertheless, previous to any perturbative decay, the evolution of the system might be dom- 
inated by non-perturbative effects as those of "preheating", i.e., the parametric amplification 
of the fluctuations in a background of oscillating fields [27,28] . Through parametric resonance, 
induced by the time dependent effective mass term in the evolution equations, field mode ampli- 
tudes can grow exponentially in time within certain resonance bands in A:-space. The paramet- 
ric resonance can be present whenever there is a non-adiabatic change in the effective masses, 
\mi{t) / mi{t)\ > 1. Thus, preheating become a more effective mechanism of transferring energy 
than the standard perturbative decay of the fields, in this case from the oscillating background 
fields to their quantum fiuctuations. The question is whether it also provides a different and 
efficient non-adiabatic source for the curvature perturbation, with super-Hubble {k <^ aH) 
fiuctuations exponentially amplified during preheating [29-33]. This would translate into an 
exponential increase in the curvature perturbation. In hybrid infiation the effect can be stronger 
first due to the presence of negative effective squared masses at the end of infiation [14,34], which 
are per se a source of instabilities in the evolution equations for the fluctuations. In addition 
the "mixing" through the potential between the flelds can also enhance the resonance [14,35] 

Whether or not super-Hubble perturbations can be parametrically amplified during preheat- 
ing is a model-dependent question [32,33] that has to be studied case by case. We have then 
numerically integrated the evolution equations for fields and fiuctuations during the oscillatory 
period, in order to see the effect on the curvature. As already mentioned, the presence of tachy- 
onic instabilities in the effective masses, combined with the fact that the fields are oscillating 
with large amplitudes and the effect of the expansion is negligible, could give rise to a strong 
parametric resonance effect just in a few oscillations. This can be seen in Fig. (3), where we have 
let the fields oscillate. It can be noticed that the resonance for the super-Hubble fiuctuations 
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Figure 3: Shown is the amplitude of the spectrum of the field fluctuations versus Nose = Ai(K(/)o)/(27r), without 
taking into account the decay of the Higgs fields. The numerical integration starts when the fluctuations left the 
Hubble horizon, i.e., k = H^^a^, roughly 60 e-folds before the end of inflation. We have taken the values of the 
parameters as before: 4>o = 10^^ GeV, = 2(^o, k0o = 1 TeV, Ax = 1.5-\/2k(/)o, = ruh = O.IH and /i* = 247 
GeV. 



only sets in once the Higgs start oscillating with an amplitude comparable to that of the others, 
disturbing the straightline trajectory followed until then by N and 0. Within that trajectory, 
their large scale field fluctuations with k <^ aH follow the same evolution equations as their 
background fields, without amplification [17]. 

Although it may seem that the amplitude of the curvature perturbation will increase beyond 
control, the above example was not realistic as we still need to take into account the effect of 
the Higgs "decay" . We remark that we used the term "decay" in an ample sense, meaning the 
transfer or dissipation of energy into radiation. In general, parametric resonance occurs within 
certain fc-bands, and if not other, the backreaction of the quantum fluctuations will tend to shut 
off the resonance. Being a non-perturbative effect, it is very difficult to study the the whole 
process without resorting to numerical lattice calculations [34,36], which are far beyond the 
scope of this paper. Analytical or semi-analytical calculations in k space have been carried out 
using the Hartree-Fock or large N approximations [37] . Using those approximations, the overall 
effect of the quantum fluctuations would be analogous to introducing "friction" in the evolution 
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equations [38,39]. Therefore, we model this effect by using a constant parameter F^, that is, a 
friction-hke term in the Higgs evolution equation [40,28], 

h + 3Hh + Vh^ -Thh . (52) 

Fh gives then the typical time interval needed for the conversion into radiation to take place. The 
assumption behind is that the quantum fluctuations are first produced by the Higgs, this being 
the field with the larger coupling to massless degrees of freedom. Although we are still working 
explicitly with the background "fields" and their perturbation variables, radiation is described 
as a perfect fiuid component, with background energy density pu and pressure pr — Pr/3. Prom 
the conservation of the total (fields plus radiation) stress-energy tensor, it follows 

PR + 3H{pR + Pr) = Vhh? = Vh{ph + Ph) , (53) 

where the RHS is the source term from the Higgs. 
Similarly for the Higgs fiuctuations we have [40]: 

^2 

Qh + SHQh +{— + Vhh)Qh + Vh4>Q^ + VhnQn = -ThQh ; (54) 
a 

here for simplicity we have neglected the sub-dominant metric contributions. In addition, we 
need to introduce their perturbations 5pR and 5Pr, in particular their corresponding gauge 
invariant variables. Gauge invariant perturbations can be defined in different ways [40], and 
a summary with our choice of gauge-invariant variables and their evolution equations is given 
in Appendix A. In order to compute the amplitude of the total curvature perturbation TZ we 
directly integrate the evolution equations for the scalar field fiuctuations, and the contribution 
of the radiation to the curvature perturbation, TZr (see Eq. (80) in Appendix A), such that 

^ = E ^^^^a = { E H4>M. + {pR + Pi^)7^i^ I ■ (55) 

By changing the patter of oscillations of the scalar fields, the parameter will control 
not only what fraction of the total energy density is transferred into radiation but how much 
"preheating" is allowed. As energy is transferred from the oscillating scalar fields to radiation 
and the amphtude of the oscillations decrease, this has the effect of shifting the resonance from 

the "broad" (large growth index for the ffuctuations) regime towards the narrow (smaller growth) 
regime, until it finally dissapears. From Fig. (3) we can read that preheating built-in on a time 
interval of the order of Atpre/i('«0o) — O(IOO). Then, depending on the value of Th/Stpreh we 
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Figure 4: Same than in Fig. (3), but including the decay of the Higgs, for two values of Th. We also 
include for comparison the analogous of the amplitude of the spectrum for the radiation, define as Pr — {{pr + 

divide the parameter range into 3 different regimes, "large" , "medium" and "small" , in order to 
study the inffuence of the friction term in the final value of the curvature: 

(a) Large: ThAtpreh > 1 (r^, ^ O(k0o) ~ 0(mi)). 

The friction is too large to allow the Higgs to oscillate (the Higgs motion is overdamped). 
The motion of the other two singlets is practically undisturbed by the presence of the Higgs, and 
they oscillate in phase along the straight line trajectory Eq. (26). Thus, there is no parametric 
resonance on the large scale perturbations. In addition, only the energy density held by the Higgs 
during inflation would be transfer into radiation, negligible compared to the total energy density. 
Inflation ends with the vacuum energy being converted into that of "matter" (oscillating fields), 
and there is no further change in the curvature perturbation. The final value of the amplitude 
of the primordial spectrum is then given by the value at the end of inflation, Eq. (48). 

(b) Medium: r^Atp^eh ~ 0(10 - 1). 

The friction is large enough to quickly shift the resonance to the narrow regime, but allowing 
first the Higgs to oscillate with an amplitude similar to the others. In this range of values, the 
energy density in radiation is comparable to the one left in the singlets system. As an example, 
on the RHS in Fig. (4) where we have plotted the amplitude for the spectrum of the fields 
fluctuations for = 0.1 (k0o)- The Higgs decays in a time interval comparable to that of the 
frequency of the oscillations, and comparable to the time needed to built the resonance for the 
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field fiuctuations, so that indeed the latter hardly starts for the fields; Higgs fiuctuations are 
just converted into radiation perturbations. 

Therefore, very quickly we are left with a mixture of radiation and two oscillating singlets. 
We may think, quite correctly, that in any case given that the individual curvature perturbations 
were more or less the same at the end of inflation, and the flrst stage of the oscillations does not 
seem either to create large differences among them, there is no further large entropy/isocurvature 
perturbation to be converted into curvature one, so that the evolution of the later would not be 
much altered. However, it is not only the presence of the radiation which make the cosmic fluid 
non-adiabatic and may affect the evolution of the curvature perturbation, but the fact that the 
flelds are still oscillating not in phase with large amplitudes. Taking now the system to be made 
of 3 components, radiation, N and fleld, the evolution equation for the total curvature can be 
written as^^ [40,19]: 

n ^ Hii-cD^Ssn + m^-cD^ins-nn) 

-mcl-c%)^{n,-nr,), (56) 
IitHs 

where we use the same notation than in Eq. (31); the subscript "T" means total quantities, and 
"S" refers to the combination of the two scalar flelds, with hs = + hjsi = (j)^ + and 

ns = ^7^^ + ^7^^ ; (57) 

lis lis 

SsR is the entropy perturbation between "scalars" and "radiation" , deflned in Appendix A, Eq. 
(85); and = pa/ Pa are the sound speed for each component. For the scalar flelds they are 
given in Eq. (32), and they are time-dependent, fast oscillating functions. Therefore, even if 
on average (over many oscillations) we could consider both flelds as behaving like matter, with 
(c^) = 0, strictly speaking the last term in Eq. (56) only cancels out when Vjv^— K/>iV = 0, which 
is no more than the condition for the straight-line trajectory. Pushed away from that trajectory 
by the Higgs, the last term in Eq. (56) introduces another non-adiabatic source in the equation, 
different from the one in the flrst fine due to the presence of the radiation. Both together can 
still enhance the curvature perturbation during the long reheating period that follows inflation. 
Until the radiation is redshifted away, we are left again only with "matter" , and the curvature 
remains constant thereon. 

(c) Small: ThAtpreh < 0(1). 
^^We neglect for simplicity terms that go like {k/aH)'^ <C 1. 
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There is still partial preheating, while the system continuosly shifts from broad to a narrow 
resonance, and the fluctuations are initially amplified. In addition, small values of Th are 
more efficient in "extracting" energy from the oscillating fields, so that by the time the Higgs 
oscillations disappear we end practically with a radiation dominated Universe, with pn/ px — 0.9. 
As an example, on the LHS in Fig. (4) where we have plotted the amplitude for the spectrum 
of the fields fluctuations^^ for Fft, = 6 x 10~'*(k(/)o). Still there is enough energy left for the scalar 
fields to oscillate, so the last term in Eq. (56) cannot be neglected. 

In Fig. (5) it is shown the spectrum of the total curvature perturbation and that of the 

radiation, for "medium" and "small" values of F/^, but using now a more cosmological time scale 

ln(a/ao), with set when reheating/preheating (oscillations) starts. The smallest F/j, the larger 

the initial amplification of the fiuctuations and then the curvature. After roughly half e-fold, the 

evolution follows a more regular pattern, with the curvature increasing as 7^ ~ ln(a/ao), until 

the radiation is diluted away and we recover TZ ~ Constant. For example, for F/j = 0.1 (k0o), 

we have {pr/ Pr)"^"'^ — 0.78, and then the radiation becomes subdominant by ln(a/ao) — 3.5, 

"'^^Numerically we found that for Th — 5 x 10^'*(k(/)o), and the other value parameters as given in the figures, 
the resonance remains broad for too long, given too large amplitude fluctuations. 
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and the curvature perturbation levels off as can be seen in the plot. The spread in the value of 
is related to the oscillations: during those the kinetic energy of the fields passes close to 
cero, at which points the curvature is bounded by the contribution in radiation. 

To summarise, we enter the reheating phase after inflation with a non-adiabatic mixture of 
oscillating fields and radiation, which will influence the evolution of the curvature perturbation 
(see Eq. (56)). On average, it will increase among one and two orders of magnitude before the 
radiation is redshifted away and becomes constant. In addition, for small values of there is 
further enhancement of the perturbations due to initial preheating effects, which are difficult to 
quantify analytically. The final spectrum would be given by that at the end of inflation, Eq. 
(48), but multiply by a function of Th, 



which depending on the value of Th can be as large as lO'""', and it could then compensate for 
the suppression of the initial spectrum of N and h during the phase transition at the end of 
inflation. 

Before ending this section, some remarks about some of the assumptions involved in the 
final value of the curvature perturbation. First, as already remarked in section 2.3, the value 
of the Higgs vev is such that it never interferes with the evolution of the other fields during 
inflation, see Eq. (24). Its evolution plays an important role during the preheating stage, but 
not before. Were this not the case, and the value of the Higgs background field h is such that 
the combination Xh became larger than H before the critical point this would signal the end of 
inflation, at which point wc would have already TZ ~ TZn — T^h — T^tp, with a different (smaller) 
suppression factor than that given in Eq. (47). In this case h would be the field with the 
increasing kinetic energy, so that it would be directly the isocurvature perturbation between the 
inflaton and h which sources the total curvature perturbation. Nevertheless, because in this case 
the initial value of the Higgs vev would be larger, its initial curvature value at horizon crossing 
becomes smaller, and we do not obtain a larger amplitude for the curvature at the end of the 
phase transition. 

For the range of values /i* for which Eq. (48) holds, the final value of the spectrum will 
scale as the inverse of the Higgs' vev, or almost. For "medium" values of Th this is the case, 
as reducing or increasing the value of h^, has null or little effect on the patter of oscillations. 
However, for "small" values of Th for which preheating is at hand, any change in the parameters 
may affect the final value, by shifting the resonance from one band to another, from broad 




(58) 
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to narrow. This will modify the initial "boost" in the perturbations, but not the afterwards 
behaviour TZ oc ln(a/ao). 

6 Discussion 

Like in the original curvaton model, the final primordial amplitude of the spectrum is con- 
trolled by the Higgs vev /i* instead of the inflaton vev at the time of horizon crossing. This 
is achieved through the field in hybrid inflation, which mediates between the inflaton and 
the Higgs /curvaton flat direction. In addition, due to dynamics of the system, the curvature 
perturbation changes first during the phase transition, which enhances the initial amplitude at 
horizon crossing by a factor ~ (^00c/''^p)^''^(0*/^*); and then through preheating effects, which 
depend on the parameter Th but nevetheless give rise to an enhancement of the spectrum. For 
small values of Th, both effect can cancel each other such that one recover a similar curvaton 
scenario result, with H^/rjJ^h*. However, even in this simplify estimation, the different 

origin of the spectrum compared to the more standard curvaton scenario can be seen in the 
factor r]J^ . That is, the amplitude of the spectrum depends to some extend on the inflationary 
dynamics, in particular on the mass of the inflaton fleld instead of the mass of the Higgs fleld. 
As far as the Higgs evolution does not backreact onto the others until the flelds start oscillating 
(see Eq. (24)), it is only its vev at horizon crossing which sets the scale for the amplitude. 

On the other hand, the evolution of the Higgs fleld and not the others will control the spectral 
index ^5 of the primordial spectrum. 

Although in the flnal value of the spectrum we have to take into account the function F[r/i]re/i 
which numerically may depend on the value of /i*, this does not bring any additional k de- 
pendence into the spectrum. That function parametrise the effects on the large scale spectrum 
during the oscillatory phase of the flelds, which depends largely on the patter of the oscillations 
and modify equally all the wavelenths with k <^ aH by the time of reheating. Given than in 
this model the variation of the Hubble parameter during inflation is practically negligible, with 
2H /H"^ ~ r]'^{4>c/ nT-py : and using dln/c ~ Hdt, we simply have: 

ns-lc^2rjhc^^, (60) 

and practically no running of the spectral index. The recent WMAP data [41], combined with 
other CMB experiments, prefers a red tilted spectrum with ns — 0.96 ± 0.02. This leads to 
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the constraint \rih\ < 0.03, which would only require ruh slightly smaller than i?, ~ 0.3if. 
Strictly speaking, the spectrum is red tilted when r]h < 0, i.e., m\ < during inflation. Although 
up to now we have been implicitely assuming the squared Higgs mass parameter in the potential 
to be positive, none of the results regarding the amplitude of the spectrum would change had 
we taken initially the opposite sign. 

7 Summary 

We have presented a variant of the curvaton scenario in SUSY hybrid inflation where: (a) it 
is coupled to the infiaton through the mediating field in hybrid inflation; (b) instead of being 
a late-decaying scalar, the curvaton decays before the inflaton, already at the beginning of the 
reheating period. The essence of the scenario is the same: isocurvature (relative) perturbations 
between the fields originated during infiation can be converted into adiabatic ones later on. 
Thus the primordial spectrum of perturbations is given in terms of the parameters related to 
the curvaton field instead of those of the infiaton. The conversion occurs when the curvaton 
energy density becomes dominant or comparable to that of the infiaton. 

However, in our scenario it is the fact that the fields are coupled which first induced the 
conversion of the isocurvature perturbation between the inflaton and the N fleld into the total 
one, previous to any decay. First, the coupling to the "curvaton" gives rise to an curvature 
perturbation given at the time of horizon crossing by the value of the curvaton one. During the 
phase transition from the false vacuum towards the global minimum, this contribution comes 
to dominate the total curvature when the kinetic energy of the N field grows until becoming 
comparable to that of the infiaton (see Fig. (2)). Therefore, the initial amphtude of the total 
curvature perturbation at the time of horizon crossing, ~ H^/rj^pcf)^, gets enhanced by a 
factor of the order of (77</,0c/"^p)^''^(0*/^*), already at the end of infiation. 

Second, the curvaton coupling to the other fields is relevant during the period of oscillations 
and it gives rise to preheating or parametric amplification of the large scale fiuctuations. In 
standard SUSY hybrid inflation, the background flelds oscillate along a practically straight-line 
trajectory. This regular behaviour prevent any ampliflcation of the fluctuations with wavenum- 
bers k <^ Ha. Including the coupled curvaton directly perturb the background trajectory, 
which immediate effect is the amplification of the large scale fiuctuations, and consequently the 
amplitude of the curvature perturbation. Without a detail numerical computation including 
backreaction and rescattering effects of the modes (large and small scale) produced, it is diffi- 
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cult to estimate how large could be this effect. However, heuristically we expect the curvaton to 
quickly decay in a few oscillations, and this would ensure that the resonance ends by the time 
the curvaton dissapears, if not before. We model then this behaviour by introducing a "decay 
rate" term in the equation of motion of both curvaton background field and fluctuations, which 
simply parametrises the energy transfer rate between the curvaton and radiation. In a sense 
this is equivalent to take into account the backreaction effect of the fluctuations. Friction or 
dissipation effects are ultimately due to quantum fluctuations, which modify the propagation of 
the fields. The approximation behind our equations is that the dominant effects are due to the 
fluctuations generated by the curvaton, with the largest decay rate, but not directly coupled to 
the singlets, and then we do not consider any other direct backreaction effect in the evolution 
equations for the singlets. Within this approximation, the smaller the decay rate parameter is, 
the larger the initial preheating and the larger the final value of the amplitude of the curvature. 
Therefore, depending on the value of F^, the amplitude enhancement during preheating can be 
larger than during the phase transition. 

This general behaviour of the curvature perturbation docs not depend on model parameters, 
but on the fact that we have coupled fields. However, the model introduced in Section 2 was 
motivated by particle physics considerations, which lead us to identify the "curvaton" D-flat 
direction with the Higgs of the MSSM. The N field then originates the term in the super- 
potential, needed for electroweak symmetry breaking. We also impose an U{1) Peccei-Quinn 
symmetry wich solves the strong CP problem. Those physical requirements fixed uniquely the 
values of the parameters in the model, as given in Section 2. The numerical results in the Fig- 
ures have been obtained using those values, mainly k0o ~ — 1 TeV and 0o — (t>c — 10^^ 
GeV. Although the vev of the Higgs D-flat direction during inflation can be considered a free 
parameter, it looks natural to assume a value not that far from its value at the electroweak sym- 
metry breaking minimum, i.e., of the order of 1 TeV. For that choice, it is clear that the value 
of the spectrum at the end of inflation, Eq. (47), although larger than at the time of horizon 
crossing, it is still far below the observational bound P^"^ ~ 5 x 10~^. However, as we have 
seen, preheating effects can further enhanced this value. Although numerically the final value 
will depend on the particular value of the parameter F/j, i.e., on how fast energy is transferred 
to the radiation, it can be seen in Fig. (5) that for example for F/j ~ 1 GeV the amplitude is of 
the correct order of magnitude. 

In conclusion we have discussed the potentially important role played by preheating in certain 
variants of the curvaton mechanism in which isocurvature perturbations of a D-flat (and F-flat) 
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direction become converted to curvature perturbations during reheating. We have analysed 
the transition from inflation to reheating in some detail, including the dynamics of the coupled 
curvaton and inflation flelds during this transition. We have discovered that preheating could be 
an important source of adiabaticity where parametric resonance of the isocurvature components 
amplifies the super- horizon fiuctuations by a significant amount. As an example of these effects 
we developed a particle physics motivated model which we recently introduced in which the 
D-flat direction is identified with the usual Higgs field. Our new results show that it is possible 
to achieve the correct curvature perturbations for initial values of the curvaton fields of order 
the weak scale. In this model we have showed that the prediction for the spectral index of 
the final curvature perturbation only depends on the mass of the curvaton during inflation, 
where consistency with current observational data requires the ratio of this mass to the Hubble 
constant to be < 0.3. 
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Appendix A 

In this Appendix we summarise our conventions and notation for the perturbations and their 
evolution equations. We work with gauge invariant (GI) quantities, to first order in linear 
perturbation theory. 

• Definitions: Linear scalar perturbations of the metric are given by the line element: 

ds^ = -(1 + 2A)dt^ + 2aBidx'dt + a^[{l - 2ip)5ij + 2E^ij]dx'dx^ , (61) 

where Bi = iki/kB, and E^ij = {—kikj/k'^)E for scalar perturbations; ip is the gauge-dependent 
curvature pertubation. On the other hand, the GI variable 

$ = --0 + Ha{B - aE) , (62) 

correspond to the curvature perturbation in the longitudinal (or zero-shear) gauge. 
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The perturbations of the stress-energy tensor are: 

ST^ = -5p, (63) 

5T^ = (p + P)(vj-Bj), (64) 

STi = -ip + Py, (65) 

ST] = SPS] + PIVj, (66) 

where Iljj is the anisotropic stress taken to be zero hereon. The comoving curvature perturbation 
is then given by: 

n = ^+^{v-B) = -^ + ^V, (67) 
k k 

where V — v — akE is the GI velocity of the fluid. 

Quantum scalar fields are decomposed as ^ai^-i t) — 4>a{t)'^^4>a{xi t), where 5(f)a{x, t) are the 
quantum fluctuations of the fleld around background values (f>a- For a multi-scalar fleld system, 
we have then: 

{p + P){v-B) = -J2^(3ScPf,, (68) 
a ^ 

and the comoving curvature perturbation: 

p-h-r^ p-\- 1- ^ 

where Qp are the Sasaki-Mukhanov GI variables: 

Qa = S(Pa + (70) 

• Evolution equations: Let us consider a multi-scalar system, with the fields interacting 
through the potential V{(j)a)- The derivaties of the potential with respect to the fields are 
denoted by Va — dV/d(f)a, and Va/s — d'^V/d(f)ad(l)p. We also allowed source terms Sa in the 
evolution equation for the background fields, 

^a + 3H^^ + Va^-Sa, (71) 

subject to the constraint Y.a Sa4>a = 0, as required by the conservation of the total stress-energy 
tensor. The corresponding GI source perturbation is defined as: 

VSa = SSa+^^, (72) 
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and the evolution of the fields perturbations is given by [40] : 

^2 

Qa. + SiJQa + ^(3q+ V'a^Q/? 
a ^ 



(73) 



where 



>l = yl + 



1 



2 p 



7^. 



In the absence of source terms, the above equation reduces to the more familiar one [25]: 



H ■ 



Q/3 



(74) 

(75) 
(76) 



The other particular case we have considered in this paper is when only one of the fields 7 
decay into radiation pr, such that 

'<i)^ + ?,H^^ + v^ = -r^0T,, (77) 
pr + ?,h{pr + Pr) = r^02 = r^(p^ + p^) . (78) 

The GI radiation perturbations associated to the energy density, prAqr, and comoving curva- 
ture, TZr, are then: 



Pr\r = opR + PR— , 

^ Ha,, 
TZr = -$ + — Vr. 



(79) 
(80) 

The variables Qa and Pa^ga a-re the fiuctuations of the field and energy density respectively on 
uniform curvature slices of space-time. For a scalar field, we have: 

(81) 



Pa^ga = "^a^ + 4>aQa + V^Qa ■ 

The evolution equations for the fiuctuations can be written as: 

.. . y^2 

a 

+ (20Q + Fq^q)— - 2-^0a— + <j>a(c% - i)PRAgR + ^a^j'j^jQ'y^ , 



3H I X ^ 



(20a + 



(82) 
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where we use the shorthand notation ha — Pa + Pa- The equations for the radiation variables 
are given by: 

H. 

h' 



{hRTZR)- = -—hR{n-nR) + HclpRAgR-3HhRnR + rJ^^n^ (83) 



k 

(pRAgR)- = -3H{l + cl)pRAgR--hRVR 

a 

+r^0^(2g^ + -^0^7^) . (84) 

The entropy perturbation Sap introduced in Eq. (56) is defined as the difference between 
the energy density perturbations: 

Pa^ga P(3\(3 

^a/3 = —r 7 , (85j 

tla tip 

which for non-interacting components is equivalent to 

= -3H f ^ - ^) . (86) 

\Pa PP J 

Appendix B 

Following Ref. [17], for a multi-scalar system the field perturbations can be rotated in field 
space into the "adiabatic" fiuctuations along the background trajectory, which originates the 
total curvature perturbation, and the "entropy" components, which are ortogonal to the fields 

trajectory. Suppose we have numbers of fields an let us denote by a the adiabatic field and 
6i the — 1 angles that parametrize the trajectory, will be the tangent vector to the field 
trajectory and vg^ the N — 1 ortogonal vectors. Denoting by Q the A^-component vector with 
the fiuctuations of the original fields, the adiabatic and the — 1 entropy fiuctuations are given 

by 

Qa^Va-Q, QOi^VQi-Q, (87) 

with 

n = , (88) 

(7 

where (7^ = J2a=i,N ^a- The evolution of the curvature for scales k <^ aH can be written as 

7t~— Yl ^iQoi, (89) 

i=l,iV-l 

This way, it is clear that the evolution of the the total curvature depends only on the entropy 
fiuctuation fields, and the curvature of the background trajectory in field space. As remarked 
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in [17], if the fields move in a straigfit line in field space, 9i — 0, the curvature perturbation on 
super- Hubble scales remains constant 

In our scenario, with 3 scalar fields, the change in the curvature during the phase transition 
can be related to the change in the trajectory in the <p — N subspace, such that the corresponding 
9i becomes non-negligible and dominant. As far as the Higgs does not contribute, the infiaton 
and N will end moving in a (practically) straight line, the system effectively reduces to that 
of one scalar field, and IZ ~ Constant. Later on, when the Higgs "curves" the trajectory, the 
evolution of TZ is modified again. 
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